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Summary 


Although  normal  stress  differences  in  liquids  have  conventionally  been  associated  with 
polymers,  aspects  of  rheological  behavior  in  lubricated  concentrated  contacts  suggest  that  normal 
stress  difference  may  be  significant  in  even  low  molecular  weight  liquids  sheared  under  high 
pressure  and  high  shear  stress.  A  torsional  flow  rheogoniometer  was  constructed  for  use  at  high 
(300  MPa)  pressure.  Four  typical  liquid  lubricants  were  investigated,  including  one 
polymer/mineral  oil  solution.  Shear  stress  and  Ni  -  N2  are  reported  as  functions  of  shear  rate. 
The  effect  of  pressure  variation  is  reported  for  two  liquids.  Results  are  compared  with  predictive 
techniques  and  a  molecular  dynamics  simulation.  Simple  low  molecular  weight  lubricant  base 
oils  can  generate  measurable  and  significant  normal  stress  differences  when  sheared  at  high  shear 
stress. 
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Introduction 


Viscoelasticity  in  liquids  has  conventionally  been  associated  with  high  polymers.  The 
related  phenomena  of  shear-thinning  and  normal  stress  differences  in  simple  shear  are  not 
ordinarily  associated  with  simple  low-molecular  weight  liquids.  The  exception  would  be  the 
lubrication  of  concentrated  contacts  with  typical  liquid  lubricants.  The  assumption  of  Newtonian 
behavior  for  the  lubricant  film  nearly  always  leads  to  large  errors  in  the  prediction  of  traction 
(friction).  Early  models  of  elastohydrodynamic  lubricant  rheology  invoked  viscoelasticity 
(Dyson,  1970  and  Chow  and  Saibel,  1971),  the  assumption  being  that  the  high-pressure  of  the 
film  sufficiently  raised  the  glass  transition  temperature  of  the  liquid  to  make  elasticity  important. 
The  shear-thinning  aspect  brought  theory  into  agreement  with  experiment.  Shear-thinning 
models  have  been  refined  to  the  extent  that  for  some  cases,  traction  curves  can  be  accurately 
calculated  from  flow  curves  obtained  with  rheometers  at  high-pressure  (Bair,  1994).  The  second 
consequence  of  viscoelastic  response,  the  normal  stress  difference,  has  been  neglected  in  studies 
of  concentrated  contact  lubrication  -  possibly  because  Tanner  (1967)  argued  that  to  affect  the 
load  capacity  of  a  lubricant  film,  the  ratio  of  the  first  normal  stress  difference  to  the  shear  stress 
must  be  significant  when  compared  with  the  length  to  thickness  ratio  of  the  film.  His  argument 
was  not  fully  accepted  at  the  time. 

Recent  developments  have  motivated  a  fresh  look  at  normal  stresses  in  liquids  under 
pressure.  Polymer  solutions  in  concentrated  contact  have  been  found  to  build  films  significantly 
thicker  than  would  be  expected  from  the  shear  viscosity  of  the  solution  when  the  films  are  very 
thin  (Cann  and  Spikes,  1994).  The  first  normal  stress  difference  acts  to  augment  the  load 
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carrying  capacity  of  a  liquid  in  shear  and  polymer  solutions  are  well  known  to  generate  a  large 
first  normal  stress  difference.  Normal  stresses  have  been  obtained  for  molecularly  thin  liquid 
films  sheared  between  atomistically  smooth  boundaries  (Reiten,  et  al.,  1994).  A  non-equilibrium 
moleeular  dynamics  simulation  of  a  liquid  lubricant  sheared  at  high-pressure  has  predicted 
(Berker,  et  al.,  1992)  a  first  normal  stress  difference  which  exceeds  the  shear  stress.  The 
orientation  of  shear  bands  which  are  observed  in  liquids  under  pressure  and  high  shear  stress 
(Bair,  et  al.,  1993)  can  so  far  only  be  explained  by  a  first  normal  stress  difference  which  is 
comparable  to  the  shear  stress  (Bair  and  Winer,  1993a).  Normal  stress  differences  may  be 
necessary  to  provide  the  geometric  instability  which  precipitates  slip  along  shear  bands  (Lee,  et 
al.,  1994)  in  liquid  lubricants. 

Background 

The  experimental  flows,  which  are  used  to  measure  the  viscometric  functions  of  liquids, 
reduce  at  least  at  some  small  scale  to  a  steady  simple  shearing  which  is  shown  in  Figure  1 .  Flow 
is  in  direction  1  with  the  velocity  gradient  in  direction  2.  The  shear  rate,  y  ,  is  the  magnitude  of 
the  velocity  gradient  and  is  twice  the  non-zero  component  of  the  deformation  rate  tensor,  d|2. 
From  symmetry  the  only  non-zero  shear  stress  is  t  =  a,,  =  Oj,  and  must  be  an  odd  function  of 
Y .  The  mean  mechanical  pressure  is  p  =  -(a,,  +  +  <^33)/ 3  and  is  not  set  by  these  flow 

kinematics.  Only  the  differences  between  normal  stresses  can  be  determined  from  the  flow  and  it 
is  known  from  symmetry  that  they  must  be  even  functions  of  shear  rate.  The  first  normal  stress 
difference,  Ni,  is  as  shown  in  Figure  1,  the  difference  between  the  normal  stresses  in  the  flow 
direction  and  in  the  velocity  gradient  direction.  The  second  normal  stress  difference  is  the 
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difference  between  the  normal  stresses  in  the  velocity  gradient  direction  and  the  neutral  direction. 
We  have 


T  =  0,2  =  ^1(7)7 

(1) 

N,  =  a„-022  =  vt/,(Y)y^ 

(2) 

N2  =cT22-a33  =  11/2(7)7' 

(3) 

where  q  is  the  coefficient  of  apparent  (or  effective)  viscosity  and  v|/i  and  v|/2 

are  the  first  and 

second  normal  stress  coefficients  respectively. 

The  large  body  of  experience  obtained  with  polymeric  liquids  tells 

us  that  all  three 

coefficients  (the  viscometric  functions)  are  monotonically  decreasing  in  7  and 

0 

11 

0 

11 

(4) 

Vi(7-^0)= 

(5) 

¥2(7^0)  =  V|/ 20 

(6) 

in  the  terminal  regime.  In  general  Ni  is  positive,  N2  is  negative  and  |N,|«|N||  .  We  will  make 

use  of  the  latter  at  times  to  justify  use  of  the  Weissenberg  hypothesis  that  N2  =  0. 

The  constitutive  equation  for  a  Newtonian  liquid  is 

=  -  Sy  Pt  +  2  p  d,.  +  p,  5y  d,k  (7) 

where  pi  is  the  thermodynamic  pressure  and  ps  is  the  second  coefficient  of  viscosity.  Then  for  a 
Newtonian  liquid  q  s  p,  Ni  s  0  and  N2  =  0.  Notice  that  the  thermodynamic  pressure  which  is 
obtained  from  an  equilibrium  equation  of  state  is  not  necessarily  equal  to  p.  Indeed,  for 
compressible  viscoelastic  liquids  with  normal  stress  differences,  a  precise  definition  of  pressure 
is  beyond  current  understanding  (Ko  and  Bogue,  1 987). 
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The  Mohr-Coulomb  criterion  is  useful  for  describing  the  orientation  of  shear  bands  in  a 
simple  shear  experiment  (Bair  and  Winer,  1993a).  Slip  will  occur  on  a  plane  for  which  the  shear 
stress  has  reached  a  critical  value  which  follows  a  linear  variation  with  the  compressive  normal 
stress  on  the  slip  plane.  At  high-pressure  the  critical  stress  becomes  nearly  proportional  to 
pressure  with  the  portionality  constant  being  an  internal  friction  coefficient  which  is  a  material 
property.  Three  idealizations  of  liquid  response  are  shown  in  Figure  2  along  with  the  resulting 
shear  bands.  The  horizontal  lines  are  solid  boundaries  to  a  liquid  film.  The  lower  boundary  is 
stationary  and  the  upper  boundary  moves  with  velocity,  v.  The  top  case  shown  in  Figure  2a  is  for 
a  pressure  insensitive  liquid  (internal  friction  coefficient  =  0)  with  zero  first  normal  stress 
difference,  Ni.  The  shear  bands  are  aligned  with  the  principal  shear  strain  rate  axes  as  shown. 
The  second  case  (Figure  2b)  is  for  a  pressure  sensitive  liquid  and  Ni  =  0.  The  shear  bands  are 
rotated  off  of  the  principal  shear  directions  in  a  sense  which  will  reduce  the  compressive  normal 
stress  on  the  planes  of  the  shear  bands.  In  Figure  2c  the  pressure  sensitive  liquid  has  an  extra 
tensile  stress  in  the  flow  direction,  i.e.,  Ni  >  0.  Effectively,  the  principal  shear  stress  axes  and  the 
shear  bands  have  been  rotated  from  the  principal  shear  strain  rate  axes  in  the  direction  of  material 
rotation.  Comparison  with  the  micrograph  of  shear  bands  (Figure  2d)  shows  that  the  complete 
description  requires  both  pressure  sensitivity  and  N  i  >  0. 

Instruments 

The  rheometers  used  in  this  study  were  of  two  types.  Measurements  of  shear  stress  as  a 
function  of  shear  rate  in  Couette  flow  were  performed  with  a  High-Pressure  High-Shear  Stress 
Viscometer  which  was  described  in  detail  elsewhere  (Bair  and  Winer,  1993b).  Measurements  of 
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both  shear  stress  and  normal  stress  difference  in  torsional  flow  were  obtained  with  a  new  High- 
Pressure  Rheogoniometer.  This  instrument  was  constructed  by  replacing  the  concentric  cylinder 
pair,  torque  transducer  and  closure  in  the  pressure  vessel  of  the  Couette  device  by  a 
rheogoniometer  cartridge  which  is  shown  in  Figure  3.  In  the  High-Pressure  Rheogoniometer  the 
pressurized  sample  liquid  occupies  the  entire  space  around  the  cartridge.  The  liquid  sample  is 
sheared  in  the  gap  between  parallel  planes  formed  by  the  end  of  a  rotating  spindle  shaft  above 
and  a  fixed  cap  below.  The  radius,  R,  of  the  shear  gap  is  2.39  mm.  The  spindle  shaft  is  made  of 
high-thermal  conductivity,  high-strength  W-Cu  composite  and  is  supported  by  two  ball  bearings 
with  low  viscosity  diester  lubricant  sealed  within  a  housing  shown  at  the  top  of  Figure  3.  The 
cap  was  constructed  of  W-Cu  composite  and  attached  through  a  non-rotating  ball  and  socket  joint 
to  a  two  channel  dynamometer.  This  dynamometer  provides  measurement  of  torque,  Q,  and 
thrust  force,  F,  on  the  cap.  It  consists  of  an  aluminum  cylindrical  shell  of  0. 1  mm  thickness  with 
two  full  strain  gauge  bridges  applied.  The  dynamometer  is  coupled  to  the  spindle  housing 
through  a  spring  preloaded  compound  screw  which  provides  adjustment  of  the  gap.  h.  by  25  pm 
per  screw  revolution.  The  strain  gauge  signals  are  passed  through  the  closure  plug  by  six 
electrodes.  The  counter  surfaces  of  the  shear  gap  were  lapped  against  one  another  with  diamond 
paste  until  the  contact  patch  was  uniform. 

Torsional  flow  between  rotating  parallel  planes  provides  measurements  of  viscometric 
functions  which  are  free  of  the  controversy  associated  with  some  other  techniques.  The 
calibration  relations  have  been  derived  many  times  and  can  be  found  in  most  rheology  texts. 
They  are  written  for  conditions  at  the  gap  edge  (maximum  shear  rate,  y  =  QR/h)as 
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(8) 


T  = 


3Q 

27iR^ 


^  ,  1  dlogQ 
3  dlogQ 


N, -N,= 


2F 

tiR^ 


^  ^  1  dlogF 
2  dlogQ 


(9) 


where  Q  is  rotating  rate  in  radians/s.  The  derivatives  in  the  brackets  provide  a  correction  for  the 
nonuniformity  of  shear  rate.  Note  that  the  normal  stress  difference  which  is  measured  is  Ni  -  N2. 
If  the  magnitude  of  N2  is  small  compared  to  Ni,  then  this  measurement  can  be  compared  to  N| 
without  significant  error. 

The  thrust  force,  F,  is  due  to  a  directly  acting  normal  stress  superimposed  on  a  pressure 
profile  which  increases  with  decreasing  radius  along  the  shear  gap.  The  normal  stress  differences 
produce  a  tensile  hoop  stress  which  generates  the  radial  pressure  profile.  The  pressure  which  we 
directly  measured  (with  a  gauge)  is  the  pressure  of  the  liquid  surrounding  the  shear  gap,  ps.  The 
boundary  condition  at  the  edge  of  the  gap  is  ps  =  -033.  From  the  definitions  of  mean  mechanical 
pressure  and  normal  stress  differences,  at  the  edge  p  =  ps  -(Ni  +  2N2)/3.  A  precise  knowledge  of 
experimental  pressure  would  require  separate  measurement  of  Ni  and  N2.  In  this  study  (N|  - 
N2)/3  was  never  greater  than  2%  of  ps  so  we  will  approximate  p  «  ps. 

The  dynamometer  was  designed  to  provide  similar  sensitivity  to  N)  -  N2  and  t  in  terms  of 
signal  volts  per  unit  stress.  Cross-talk  was  found  to  be  less  than  3%. 


Experiment 

Four  liquids  which  have  been  used  as  concentrated  contact  lubricants  were  the  subject  of 
this  study.  Three  of  the  liquids  are  base  oils  and  these  are  some  of  the  most  thoroughly 
characterized  (Bair,  1994  and  references  therein)  liquids  in  the  tribology  literature.  They  are  143 
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AD  -  perfluoropolyalkylether,  LVI  260  -  mineral  oil,  and  5P4E  -  polyphenyl  ether.  A  polymer 
blend,  LF  5346,  was  also  included,  because  polymer  solutions  are  well  known  to  display  normal 
stress  differences  at  atmospheric  pressure  and  they  are  an  important  class  of  lubricants.  The 
polymer  is  polybutene  with  number  average  molecular  weight  of  25,000  and  weight  average 
molecular  weight  of  50,000  and  is  blended  with  80  percent  mineral  oil  with  kinematic  viscosity 
of  18  mm^/s  at  38°C. 

The  operation  of  the  Couette  flow  device  has  been  described  previously  (Bair  and  Winer, 
1993b).  For  the  Fligh-Pressure  Rheogoniometer,  the  shear  gap,  h,  was  determined  by  generating 
flow  curves  in  both  Couette  and  torsional  flow  on  the  same  material  at  the  same  pressure  and 
temperature  and  shifting  the  torsional  flow  data  by  varying  h  until  they  superimposed. 

Viscous  heating  is  reduced  by  maintaining  h  as  small  as  possible.  Experiments  in 
torsional  flow  with  h  as  small  as  2  pm  gave  poor  repeatability,  possibly  due  to  significant  relative 
changes  in  the  gap,  h,  with  normal  force  because  of  the  compliance  of  the  dynamometer.  A  gap 
of  9.7  pm  was  used  for  all  reported  results  as  a  good  compromise  between  viscous  heating  and 
repeatability.  For  this  value  of  h  and  a  typical  temperature  viscosity  coefficient  and  liquid 
thermal  conductivity,  viscous  heating  will  be  significant  for  a  dissipation  of  ty  =  10'"  W  /  m’ . 
This  corresponds  to  a  Brinkman  number  of  one.  The  Couette  flow  measurements  were 
performed  with  a  gap  of  1.2  pm  and  a  significantly  higher  viscous  heating  limit  of  about 
xy  =  2  xlO"  W/m^ 

A  t)q>ical  thrust/torque  history  is  shown  in  Figure  4.  At  start-up  the  rotational  velocity 
was  increased  linearly  with  time  until  the  test  velocity  was  reached  at  50  ms.  Note  that  the 
torque  measurement  response  is  much  faster  than  the  thrust  measurement  response.  This  is  the 
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"squeeze  damping  effect"  which  results  from  the  small  increase  of  the  gap  in  response  to  the 
thrust  force  and  the  retardation  of  flow  into  the  gap  by  the  very  large  viscosity  of  the  liquid. 
When  the  direction  of  rotation  was  reversed,  the  thrust  was  unchanged  and  the  torque  changed 
sign  as  expected. 

Results 

Flow  curves  for  the  four  liquid  lubricants  are  shown  in  Figures  5-8  at  indicated  pressure 
where  the  closed  circles  represent  shear  stress  measured  in  Couette  flow  and  all  other  data  points 
are  for  torsional  flow.  All  tests  were  at  20°C.  These  are  apparently  the  first  measurement  of 
normal  stress  in  liquids  under  high-pressure.  For  the  three  base  oils;  143  AD,  LVI  260,  and 
5P4E  the  normal  stress  difference,  Ni  -  N2  was  always  less  than  the  shear  stress,  x.  For  the 
polymer  solution  the  normal  stress  difference  was  always  greater  than  the  shear  stress  for  the 
rates  of  shear  investigated.  Note  that  for  LF  5346,  N|  -  N2  attained  a  value  of  12  MPa  which  may 
be  the  highest  experimentally  measured  normal  stress  difference  for  a  liquid. 

Viscous  heating  effects  can  clearly  be  seen  for  the  torsional  flow  (open  circles)  in  Figures 
5  and  6  where  the  shear  stress  becomes  less  than  that  measured  in  the  Couette  device  for 
xy  >  10'“W/m'. 

Discussion 

Methods  for  obtaining  engineering  estimates  of  the  first  normal  stress  difference  when  the 
viscosity  function  is  known  are  available  in  the  rheology  literature  (Dekee  and  Stastna,  1986). 
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Two  are  listed  here.  The  first  is  due  to  Abdel-Khalik,  et  al.  (1974)  and  was  derived  from  a 
rheological  equation  of  state. 


¥ 


W'tIc 


^  r  nCf)-  nCs) 

2  -2 


ds 


0  S  -  Y 


(10) 


2  <k<3 


The  second  is  due  to  Gleissle  (1982)  and  follows  from  empirical  rules. 

Jn(Qo)  y 


(11) 


2<k<3 


A  recursion  formula  for  approximation  of  (1 1)  was  also  given. 


¥,  (k(y  -  Ay))  =  ¥i  (ky)  - 

yj{y-Ay)y 


(12) 


Note  that  the  dimensionless  shifting  parameter  k  has  a  different  meaning  in  equations  (10)  and 
(11)  although  the  stated  range  is  the  same.  In  the  former,  k  shifts  \|/i  and  in  the  latter  k  shifts  y  . 

The  line  plotted  through  the  shear  stress  data  in  Figure  8  represents  the  power  law  for  a 
power  law  exponent  of  N  =  0.35. 

q  =  Ay-^  (13) 


Equation  (10)  combined  with  the  power  law  (13)  was  solved  by  Abdel-Khalik  (1974)  and  is 
shown  in  Figure  8  for  k  =  3.  Inserting  the  power  law  (13)  into  equation  (11)  results  in 

=  (14) 


which  is  plotted  as  Gleissle  in  Figure  8  for  k  =  3.  These  predictions  bracket  our  measurements  of 
Ni  -  N2  for  the  polymer  solution. 
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It  is  of  interest  to  apply  these  engineering  predictions  to  low  molecular  weight  liquids 
which  have  not  previously  been  considered  in  tests  of  these  predictions.  The  viscosity  function, 
r],  is  plotted  against  shear  rate  in  Figures  9-1 1  for  the  base  oils.  Again,  solid  circles  are  Couette 
flow.  For  shear  rates  which  are  greater  than  the  experimental  capability,  equation  (14)  was 
applied  and  for  shear  rates  at  which  q  has  been  measured  the  recursion  formula  (12)  was  used  to 
give  the  predicted  first  normal  stress  coefficient,  v|/i  plotted  in  Figures  9-11.  The  values  of  k 
required  to  obtain  the  fit  to  the  experimental  data  are  shown  in  the  figures.  Remember,  the  data 
points  plotted  for  comparison  are  (N,  -N,)/ or  v|/,  -  \]i 2- 

We  know,  however,  that  the  measurements  of  viscosity  shown  in  the  figures  do  not 
represent  constitutive  behavior  for  rates  greater  than  the  rate  at  which  the  first  shear  bands  appear 
(Bair  and  Winer,  1993).  A  portion  of  the  velocity  used  in  the  calculation  of  the  velocity  gradient 
in  these  measurements  was  contributed  by  slip  along  the  shear  bands.  Therefore,  as  a  check,  in 
Figure  10  for  the  mineral  oil,  it  was  assumed  that  the  constitutive  behavior  followed  a  power  law 
above  the  arbitrary  shear  rate  of  250  s'‘  as  shown  by  the  broken  line.  The  prediction  for  yi  is 
nearly  the  same  for  the  two  approaches  within  the  experimental  range  of  measurement  of  vi  -  \\i 2 
and  we  can  conclude  that  localization  in  Couette  flow  does  not  affect  our  comparisons  -  only  the 
extrapolations. 

For  polymers,  k  varies  from  2  to  3;  however,  for  the  low  molecular  weight  base  oils 
investigated  we  found  a  value  of  0.4  to  1  to  be  more  appropriate.  This  may  be  the  result  of  a 
difference  in  molecular  mechanism  for  normal  stresses  for  the  different  materials.  The  non¬ 
linear  stretching  of  entangled  molecules  will  not  enhance  the  extra  tensile  stress  in  the  flow 
direction  for  the  low  molecular  weight  base  oils. 


12 


The  relationship  between  Nj  (or  Na)  and  t  has  been  found  to  be  a  characteristic  of  a 
particular  polymer  (eg.  Tanner,  1973).  Han  and  Jhon  (1986)  showed  that  Ni  vs.  y  isotherms  for 
various  temperatures  could  be  shifted  to  superimpose  upon  a  master  curve  by  plotting  N  i  vs.  x. 
Tanner  (1985)  showed  that  the  relation  between  Ni  and  x  should  be  a  weak  function  of 
temperature  by  time-temperature  superposition.  We  know  that  time-temperature  shifting  can  be 
generalized  to  include  pressure  as  well  as  temperature  so  that  the  relation  between  Ni  and  x 
should  also  be  insensitive  to  pressure.  We  have  plotted  Ni  -  N2  versus  x  in  Figure  12  for  LF  5346 
and  5P4E.  Also  included  are  results  at  reduced  pressure  as  represented  by  the  solid  data  points. 
The  shear  rates  for  these  reduced  pressures  were  four  and  ten  times  greater  for  the  same  shear 
stress  for  LF  5346  and  5P4E  respectively  because  the  low  shear  viscosity,  p,  was  reduced  by  this 
ratio.  However,  the  relationship  between  Ni  -  N2  and  x  as  shown  in  the  figure  is  little  changed  by 
pressure. 

Continuum  models  predict  that  in  the  terminal  regime  (Larson,  1988) 

0 

The  Rouse  molecular  model  gives  (Larson,  1988)  for  the  modulus,  G  =  nkT  so  that 

(16) 

where  k  is  the  Boltzmann  constant,  T  is  absolute  temperature  and  n  is  the  number  of  polymer 
molecules  per  unit  volume.  Equation  (16)  is  plotted  for  the  polymer  blend,  LF  5346,  in  Figure 
12  assuming  N2  =  0.  The  molecular  volume  of  5P4E  (molecular  weight,  450)  was  used  to 
calculate  n  and  the  result  is  also  plotted.  Both  results  (the  solid  lines  in  Figure  12)  may 
adequately  represent  the  low  shear  stress  behavior  of  these  liquids.  The  value  of  G  for  5P4E,  6 
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MPa,  is  two  orders  of  magnitude  lower  than  the  limiting  high  frequency  shear  modulus,  Goo, 
which  has  been  measured  by  various  techniques  (Bair,  1994).  G  should  be  the  modulus 
associated  with  the  longest  relaxation  time. 

The  orientation  of  shear  bands  in  simple  shear  was  noted  for  5P4E  and  the  ratio  of  first 
normal  stress  to  shear  stress  was  calculated  from  Mohr-Coulomb  theory  to  be  1 .26  (Bair  and 
Winer,  1993a).  The  shear  stress  was  approximately  20  MPa.  The  broken  line  in  Figure  12  is  Ni 
=  1 .26  T  assuming  N2  =  0.  Extrapolation  of  terminal  behavior  would  reach  the  predicted  ratio 
well  in  advance  of  20  MPa  shear  stress;  however  the  experimental  measurements  show  a  trend 
which  falls  short  of  the  predicted  ratio  of  Ni  to  x.  The  calculation  of  vi/,  (y)  for  5P4E  in  Figure 

11  by  equation  (11)  for  k  =  1  gives  the  result  that  at  high-shear  stress.  If  k  =  1.12 

instead,  the  ratio  of  Ni  to  x  required  by  the  shear  band  orientation  is  obtained. 

Berker  and  coworkers  (1992)  have  performed  a  non-equilibrium  molecular  dynamic 
(NEMD)  simulation  of  a  lubricant  under  concentrated  contact  conditions.  The  model  lubricant 
was  the  hydrocarbon,  n-hexadecane,  at  high-pressure  and  high-temperature.  The  simulation 
allowed  chemical  bonds  to  stretch  and  bend.  Homogeneous  constant  volume  shear  was  imposed 
and  the  viscometric  functions  were  reported  over  a  large  range  of  shear  rate.  The  constant 
volume  constraint  caused  the  mean  mechanical  pressure  to  increase  with  shear  rate  (and  shear 
stress).  However,  we  have  shown  above  that  a  presentation  of  normal  stress  difference  versus 
shear  stress  removes  or  reduces  the  pressure  effect.  The  results  of  Berker,  et  al.  (1992)  are 
plotted  in  Figure  13.  At  the  lowest  shear  rate  (and  lowest  shear  stress)  a  large  uncertainty  was 
reported  as  shown  in  the  figure.  The  present  experimental  results  for  the  hydrocarbon,  LVI  260, 
are  plotted  in  the  same  figure  and  are  consistent  with  the  NEMD  results.  The  molecular  weight 
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of  LVI  260  is  425;  about  twice  that  of  hexadecane.  The  terminal  normal  stress  difference 
according  to  equation  (16)  should  be  proportional  to  molecular  weight  for  equal  density.  Halving 
the  measured  normal  stress  difference  for  LVI  260  would  place  the  results  in  the  center  of  the 
results  of  the  NEMD  simulation. 

Normal  stress  differences  have  not  been  measured  for  low  molecular  weight  base  oils  at 
atmospheric  pressure  by  conventional  rheogoniometers.  The  reason  for  this  is  apparent  from 
equation  (16).  If  the  smallest  measurable  Ni  is  10  Pa,  the  shear  stress,  t,  must  be  at  least  5000 
Pa.  This  shear  stress  is  difficult  to  reach  with  a  relatively  low  viscosity  (~1  Pa  s)  liquid  in  a 
conventional  rheometer  without  significant  viscous  heating.  Also,  for  low  viscosity,  inertia 
effects  dominate  over  the  normal  stress  for  conventional  measurement  techniques.  High-pressure 
is  helpful  in  that  the  viscosity  can  be  very  large. 

Conclusions 

1)  A  High-Pressure  Rheogoniometer  has  been  constructed  which  is  capable  of  measuring  shear 
stress,  X,  and  the  normal  stress  difference,  Ni  -  N2  to  pressures  of  240  MPa  and  shear  rates  of 
at  least  10“^  s"'. 

2)  Simple  low  molecular  weight  liquids  generate  measurable  and  significant  normal  stress 
differences  when  the  shear  stress  is  high  and  the  measured  values  agree  with  NEMD 
simulation. 

3)  Methods  for  making  engineering  estimates  of  first  normal  stress  difference  may  be  useful  for 
lubricant  base  oils  if  the  shifting  parameter,  k,  is  less  than  the  value  used  for  polymers. 

4)  The  relationship  between  Ni  -  N2  and  x  is  a  weak  function  of  pressure. 
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5)  Measurements  of  normal  stress  differences  are  not  entirely  inconsistent  with  shear  band 
observations;  however,  normal  stress  measurements  obtained  at  the  critical  shear  stress  for 
slip  are  needed. 
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Pressure  Insensitive,  N,  =  0 


Pressure  Sensitive,  >  0 


Figure  2.  Shear  Band  Orientation  for  Mohr  Coulomb  Theory  and 
Various  Idealizations. 


Figure  3.  High-Pressure  Rheogoniometer  Cartridge  and  Closure. 
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Figure  4.  Typical  ThrustTTorque  Histories  for  Torsional  Flow. 


Figure  5.  Flow  Curves  for  Perfluoropolyalkylether.  Solid  Circles 
from  Couette  Shear.  Others  from  Torsional  Flow. 


Figure  6.  Flow  Curves  for  Mineral  Oil.  Solid  Circles 
from  Couette  Shear.  Others  from  Torsional  Flow. 


Figure  7.  Flow  Curves  for  Polyphenyl  ether.  Solid  Circles 
from  Couette  Shear.  Others  from  Torsional  Flow. 


Figure  8.  Flow  Curves  for  Polybutene  in  Mineral  Oil.  Power  Law  Fitted 
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ire  9.  Viscosity  and  First  Normal  Stress  Coefficient  versus  Shear  Rate 
for  143  AD.  The  prediction  of  Equation  (11)  is  plotted,  vj/i  -  \j/2  plotted 
IS  +. 


Figure  10.  Viscosity  and  First  Normal  Stress  Coefllcient  versus  Shear  Rate 
for  LVI 260.  The  prediction  of  Equation  (11)  is  plotted.  -  Xj/j  plotted 
as  +. 
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Figure  13.  Comparison  of  Molecular  Dynamic  Simulation 
with  Experimental  Measurements. 
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